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Observation of two new narrow mesons in the D + s tc and D + s 7Ty final states. 
Results from BaBar, Belle and CLEO. 
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The BaBar experiment has discovered a narrow state, denoted as D*,(2317) + , near 2.32 GeV/c 2 in the D + S n° invariant mass distribution 
from inclusive e + e~ interactions at center-of-mass energies near 10.6 GeV. The same experiment has hinted the presence of a second 
new state near 2.46 GeV/c 2 in the D*7t°y mass spectrum. These discoveries have triggered CLEO and Belle experiments in their own 
search for new states coupled to D + s meson. They have confirmed the existence of the £>*j(2317) + state and established the evidence for 
the second state denoted as Z)*j(2458) + . Belle has reported evidence for both states not only in continuum e + e~ annihilations but also 
in B decays. These two new narrow resonances are characterized by masses significantly lower than predicted by potential models and 
enough to close off the most natural decay modes, thus forcing isospin-violating transitions. Their discovery has been a surprise because 
it contradicts the expectations of apparently established theory. This circumstance has led to a great effort on theoretical side as well. 



1 Introduction 

Potential models of mesons made up by an heavy quark 
and a light one have had so far reasonable success in de- 
scribing the spectroscopy of D, D s , B and B s systems [ 
Q). The spectroscopy of cs states is straightforward in the 
limit of large charm quark mass [|2): the total angular mo- 
mentum of the light quark, j-l + s, obtained by summing 
its orbital and spin angular momenta, is conserved. The 
P-wave states, all of which have positive parity, then have 
7 = 3/2 or j = 1/2; after combination with the heavy quark 
spin (J = j + s/,), the former gives total angular momentum 
J = 2,1 whereas the latter gives J = 1,0. The 2 + and one 
1 + state, forming the doublet associated to 7 = 3/2, are ex- 
pected to have a small width because their dominant (OZI- 
and isospin-favoured) decays to D*K and DK respectively 
would proceed via D-wave. On the contrary, the + and the 
other 1 + state, associated to j = 1 /2, are expected by most 
potential models to have comparable masses (in the range 
2.4-^2.6 GeV/c 2 ) and to decay by kaon emission to the 
same final states, but with large widths (at least 200 -f 300 
MeV) since these decays would proceed via S -wave. 

The experimental information on the S -wave ground state 
is well established. The D* + meson is identified with the 
excited 1" triplet state at 21 12 MeV/c 2 , decaying to the D + s 
meson identified with the ground 0~ singlet state at 1969 
MeV/c 2 , mainly by photon emission. An isospin-violating 
pion emission also occurs with a rate of about 6%. The 
experimental information concerning the P-wave states is 
limited instead. The observations of the D* y (2573) + (with 
J p consistent with 2 + ) and the J p = 1 + D. s i(2536) + identi- 
fied with the j = 3/2 doublet were feasible because of their 
being narrow. However there are no experimental candida- 



*Talk at the International Workshop on QCD: QCD@Work 2003 - Con- 
versano (Italy) 14-18 June 2003 (eConf C030614). 



tes for the j = 1 /2 doublet so far; the observation of these 
two states should be complicated by their large predicted 
widths. 

The experimental results discussed in this review paper 
concern the recent observation of two states consistent with 
their identification with the (0 + , 1 + ) missing doublet; they 
are below DK and D*K thresholds and decay to D + S n° and 
D* + 7t° final states respectively. 



2 Observation by BaBar of the narrow state 

D;/2317) + 

The BaBar analysis, described in detail elsewhere [0, in- 
vestigates the inclusively-produced D+n mass spectrum 
by combining charged particles corresponding to the de- 
cays D + s — > K + K n + [ 4 1 with 7r° candidates that (1) are re- 
constructed from a pair of photons with energy greater than 
100 MeV and constrained to emanate from the intersection 
of the K + K'n + candidate trajectory and the beam enve- 
lope, and (2) are formed by performing a one-constraint fit 
to the 7r° mass. A given event may yield several acceptable 
7i° candidates, but only those that do not share any photon 
with the others are retained. To reduce backgrounds, both 
that from B meson decays and the combinatorial one from 
continuum, each K + K~tt + tt q candidate is required to have 
a momentum p* in the e + <?~ center-of-mass frame greater 
than 2.5 GeV/c . The upper histogram in Figf^a) shows 
the K + K n + invariant mass distribution for all candidates 
characterized by the peaks associated to D + and D + s mesons 
signal. 

Further background rejection, providing the lower his- 
togram of Figf^a), is achieved by the following addi- 
tional selection criteria. Quasi-two body decay modes into 
<pn + and K K + final states are selected by retaining only 
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Figure 1. The BaBar Collaboration [ 3|. (a) The K + K n + mass 
distribution for all selected candidates; additional selection crite- 
ria described in the text have been used to obtain the lower his- 
togram, (b) The mass distribution of two-photon combinations 
in the selected events, (c) The K + K~n + n° mass distribution for 
candidates in the D + s signal/sideband regions of (a) (top/shaded 
histogram respectively); the n° candidates are mass-constrained, 
(d) The D\yy mass distribution for candidates in the D + s signal re- 
gion of (a) coupled with photon pairs from the n° signal/sideband 
regions of (b) (top/shaded histogram respectively); both top and 
shaded histograms present a broad peak due to random combina- 
tions D * + y where D* + — 

those candidates with K + K~ mass within 10 MeV/c 2 of 
the 0(1020) mass or with K ti + mass within 50 MeV/c 
of the K (892) mass. Given the D + s Dalitz plot structure, 
these candidates belong to two disjoint sub-samples. More- 
over the decay products of the vector mesons 0(1020) and 
— *o 9 

K (892) exhibit the cos 6/, behaviour expected because of 
angular momentum conservation, where is the helicity 
angle; the signal-to-background ratio is improved by re- 
quiring | co sOh | > 0.5. 

In Fig^a) the signal and sideband regions are shaded 
and the D + s signal peak contains about 80, 000 candidates 
and is centered at a mass of (1967.20+0.03) MeV/c 2 (sta- 
tistical error only). Fig^b) shows the mass distribution for 
all two-photon combinations associated with the selected 
events; the n° signal and sideband regions are shaded. Can- 
didates in the D + s signal region of Fig^Ja) are combined 
with the mass-constrained n° candidates to provide the top 
histogram of Figtnc). The shaded histogram in Fig^c) 
represents combinations of candidates from the D + s side- 
band regions with the same n° candidates. In Fig^d) the 
mass distributions result from the combination of candida- 



Figure 2. The BaBar Collaboration [ 31). The D*n° mass dis- 
tribution for (a) the decay £>+ — > K + K~n + and (b) the decay 
£>+ K + K-jT + n°. Fit function and results for the D* y (2317) + sig- 
nal are described in the text. The very narrow peak near threshold 
is associated to the isospin- violating D* + — > Z)+ 7T° decay. 

tes from the D + s signal region with the photon pairs from 
the n° signal and sideband regions of Fig[TJb); the mass re- 
solution is poorer than in Fig^c) because there is no mass 
constraint for the photon pairs. In order to improve the re- 
solution the nominal D + s mass has been used to calculate 
the D + s energy. 

The top histograms of Fig^c) and Figlfld) (associated to 
signal-signal combinations) are characterized by a clear 
narrow signal at a mass close to 2.32 GeV/c 2 . In the 
corresponding shaded histograms (associated to signal- 
sidebands combinations) this signal is absent. These facts 
indicate unambigously that the peak is associated with the 
D^n° system. Moreover the signal is observed in both 
the quasi-two body D + s decay modes. Furthermore Monte 
Carlo simulations have been used to investigate the pos- 
sibility that the D* y (2317) + signal could be due to a re- 
flection from other known charmed states. No hint of any 
peak is found in the signal region. In addition no peak has 
been produced by exchanging the kaon and pion identities 
to simulate particle mis-identification effects. 

Fig|5]shows the D^n mass distribution for p*{D + s n a ) > 3.5 
GeV/c for two different decay modes: (a) — > 
K + K~n + and (b) D+ -> K + R-n + n {) . This second addi- 
tional sample of decays is obtained by adding to each 
K + K~n + combination fitted to a common vertex the 7r° can- 
didates constrained to emanate from this vertex; each re- 
sulting D + s candidate is combined with a second tt° can- 
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Figure 3. The BaBar Collaboration [ 3|. The mass distribution 
for (a) D + S y, (b) D+yy and (c) D*n°y. Details on combination 
criteria are given in the text. The shaded histograms of (b) and (c) 
correspond to D+y masses falling in the D* + signal region. The 
vertical line indicates the £)* J (2317) + mass value. 
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Figure 4. The BaBar Collaboration. (Right) Simulated data 
e + e~ —> cc including D* sJ (22>il) + . The scatter plot shows how a 
peaking background, at about 2.46 GeV/c 2 in the D+7r°y mass 
projection, can be generated by the kinematical overlap of the 
D* + -> D* s y signal band with the D* y (23 17) + -> D+tt° signal band 
obtained adding respectively an unrelated rc° and an unrelated y. 
(Left) Real data are characterized, at least partially, by the same 
kinematical cross-through. 



didate with an harder lab momentum (greater than 300 

— *o 

MeV/c). Moreover the K , K , cf> or p mass regions 
are selected for the relevant two-body subsystems. The fit 
function includes a single Gaussian for the signal and a 
polynomial function for the background distribution. The 
fit yields 1267 + 53 candidates over background for a fit- 
ted mass of 2316.8 ±0.4 MeV/c 2 with a standard deviation 
of 8.6 + 0.4 MeV/c 2 (statistical errors only), for the D + s 
decay mode (a). The Gaussian fit for the D + s decay mode 
(b) is fully consistent with the previous one since it yields 
273 ± 33 candidates with a mean of 2317.6 + 1.3 MeV/c 2 
and a width of 8.8 ± 1.1 MeV (statistical errors only). The 
systematic uncertainty in the mass is conservatively esti- 
mated to be less than 3 MeV/c 2 . 

An estimate of the mass resolution for the K + K~n + n° sys- 
tem can be provided directly from the data by fitting the 
mass distribution for D + s — > K + K~n + 7T° and the measured 
width is consistent with that of D* y (2317) + signal. A si- 
milar mass resolution can be obtained from the simulation 
of the D* y (2317) + decaying to K + K~n + n° with a genera- 
ted intrinsic width close to zero. In other words the width 
observed in the data is consistent with Monte Carlo expec- 
tations for a zero intrinsic width state. Thus the observed 
width is consistent with the experimental resolution and it 
can be concluded that the intrinsic width of the D* y (23 17) + 
must be below 10 MeV. 

A search for other D* y (2317) + decays has been per- 
formed. For this search additional background suppres- 
sion is achieved by requiring that the p* of each system 



of particles investigated must be greater than 3.5 GeV/c . 
FigGJa) shows the D + S y mass distribution obtained by com- 
bining a Z)j candidate from the signal region of FigQa) 
with photons in the selected events not belonging to a yy 
combination in the signal region of Figf^b) and having an 
energy greater than 150 MeV. A clear D* + signal is visi- 
ble but there is no indication of a D* y (2317) + production. 
No signal near 2.32 GeV/c 2 is observed in the D + S yy mass 
distribution presented in Fig|3lb) as well; again any pho- 
ton belonging to a yy combination in the n° signal re- 
gion is excluded here. The shaded histogram of Figl^Jb) 
is associated to the subset of combinations for which either 
D + s y combination falls in the D* s + signal region (defined as 
2.096 <m{D + s y)< 2.128 GeV/c 2 ). A D* y (2317)+ signal is 
missing here too; this implies the absence of a D* + y decay 
mode at the present level of statistics. 

The top histogram of Fig|3lc) shows the D+7T°y mass dis- 
tribution excluding again any photon that belongs to the 
n signal region. The subset of combinations in which 
the D + S y mass falls in the D* + signal region is represented 
by the shaded histogram of Fig|5Jc). No signal near 2.32 
GeV/c 2 is observed in both distributions. However a peak 
is visible near a mass value of 2.46 GeV/c 2 and survives 
at the D* + selection. This mass corresponds to the overlap 
region of the D* + D + S y and ZT y (2317) + -> £>>° sig- 
nal bands that produces, because of the small widths of 
both the D* + and D* y (23 17) + mesons, a narrow peak in the 
D + S n°y mass distribution. It has been checked the possi- 
bility that this peak at about 2.46 GeV/c 2 would kinema- 
tically feed-down the D* y (2317) + signal, provided that it 
could eventually be a new additional narrow state decaying 
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Figure 5. The CLEO Collaboration [0. Am(D>°) mass dif- 
ference distribution. The points represent the data; the solid 
histogram is the predicted normalized spectrum from simulated 
e + e~ — > qq events. 

to D* s + 7i°. This eventual feed-down would be produced as 
follows: e + e- -» X(2460) + X„, col7 , X(2460) + -> D*+n° and 
D* + — > D + S y where the D^tt combination could emulate a 
D* y (2317) + candidate and the residual photon would ap- 
pear to be unrelated or simply lost. Such D + S n° reflected 
peak, however, would have a root-mean-square of about 
15 MeV/c 2 , namely quite twice the width of the observed 
Z)* y (2317) + signal. The second crucial check is obtained 
by a Monte Carlo unfolding method: if the apparent sig- 
nal at 2.46 GeV/c 2 were due to a state decaying entirely 
to D* + 7r°, it would produce only one-sixth of the observed 
signal at 2.32 GeV/c 2 . 

3 Confirmation of D* 7 (2317) + and observa- 
tion of D* y (2458) + by CLEO and Belle 

CLEO [ 13 and Belle [ 6| have readily confirmed the 
D* y (2317) + state seen by BaBar, using only the (fm + quasi 
two-body D + s decay mode. In terms of the mass difference 
Am(D^7T°)^m(K + K-7T + jT°) - m{K + K-n + ), CLEO reports a 
value of Am = 350.3 il.OMeV/c 2 with 165+20 candidates 
(Fig|5}, whereas Belle reports Am = 348.7 ± 0.5 MeV/c 2 
with 761 +44 candidates (Fig|6j. These values (with sta- 
tistical errors only) are in good agreement with the BaBar 
measurement of Am = 348.4 + 0.4 Me V/c 2 [[3). Using 
their Monte Carlo simulation CLEO is able to limit the na- 
tural width of the D* y (2317) + state to be less than 7 MeV 
at the 90% confidence level (C.L.). 

As just discussed, in order to establish the existence of 
the Z)* y (2317) + , BaBar has ruled out any relevant feed- 
ing down by a state of mass 2.46 GeV/c 2 eventually pro- 
duced in addition to the D* y (2317) + . On the other hand 
BaBar has been rather careful to state that the structure ob- 
served in the D + S n°y mass spectrum is another new nar- 
row state; the complexity of the overlapping kinematics 
of the D* + D + S y and ZT y (2317) + D>° decays has 
suggested a dedicated and detailed study to evaluate the 
amount of feeding-up from the D* y (2317) + . Fig@Jnght) 
illustrates, by means of a Monte Carlo study, the kine- 
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Figure 6. The Belle Collaboration [|6). Ara(D>°) mass diffe- 
rence distribution. The points represent the data; the histograms 
are associated to D* (solid) and n° (dashed) sideband regions. 

matical cross-through: the D* + — > D^y signal band, ob- 
tained by adding a random n°, crosses the band produced 
by D* y (2317) + — > D+n plus an additional random y, thus 
generating a peak at about 2.46 GeV/c 2 in the D + s n°y mass 
projection. This pattern also characterizes, at least par- 
tially, the same scatter plot for the real data in Fig0Jleft). 

CLEO [ |3 and Belle [ |6) have analyzed the D* + n Q 
mass spectrum finding evidence for the structure in the 
2.46 GeV/c 2 region and claiming it is a new state, the 
Z)* y (2458) + . The mass difference is now defined as 
Am(D* + 7r°) = m(K + K-7T + n y) - m{K + K-n + y) and its dis- 
tribution is provided by CLEO in Fig0 and by Belle in 
FigOD CLEO has addressed the problem of feed-up in 
two different ways. By an unfolding method using simu- 
lated data they estimated a rather small amount of feed-up 
(about 9%) which is consistent with the small peak in the 
D* + sidebands (FigHb)). On the other hand CLEO fits 
the signal peak after performing a D* + sidebands subtrac- 
tion. The two methods provide consistent results in terms 
of D* y (2458) + candidates yield and the mean and the width 
values of the fitted D* y (2458) + signal peak in Am(D* + w°) 
distribution. Belle's observation [ 6| of the D* y (2458) + 
state is given in Fig|8] where there is also visible a clear 
peaking in the D* + sideband regions that shows a relevant 
level of feed-up (about 30% [ 7 1). Afterwards BaBar has 
confirmed the existence of the D* y (2458) + state [|U [|5J. 
The D* y (2458) + mass measurements reported by Belle [ 6 1 
and BaBar [|9) are in agreement: 2456.5 ±1.7 MeV/c 2 
and 2458.0+ 1.4 MeV/c 2 respectively (with combined sta- 
tistical and systematic uncertainties); however the mass 
value obtained by CLEO is two standard deviations larger: 
2463.5 + 2.1 MeV/c 2 . All three experiments have used 
Monte Carlo simulations to evaluate the detector mass re- 
solution. CLEO [ 5 1 provides a 90% C.L. upper limit of 7 
MeV on the natural width of the D* y (2458) + state. 

These two new states have been observed for the first time 
by Belle also in B decays \ 1101 . Fig|5] shows the recon- 
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Figure 7. The CLEO Collaboration Am(D*+7T°) mass diffe- 
rence distributions, for (a) D* + signal region and (b) D* + sideband 
regions. The signal yield obtained by the fit is 55 ± 10 candidates. 
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Figure 8. The Belle Collaboration [0. Am(D* + 7r°) mass diffe- 
rence distributions showing the evidence for £>* y (2458) state. The 
shaded histogram is associated to the D* + sidebands. The signal 
yield obtained by a fit (not shown here) is 126 ± 25 candidates. 

structed D*n , D* + n° and D\j mass spectra for candida- 
tes whose mass and beam energy constraints are consistent 
with the decay B — > DD* sJ . The signal in Fig|9jc) represents 
the first observation of the radiative decay Z)* y (2458) + — > 
D + S y. From the two B — > Z)Z)* y (2458) branching fraction 
measurements, Belle determines the ratio 



S(D; y (2458)^D iy ) 
S(D* y (2458)^D*7r°) 



0.38 ±0.1 l(stat.) + 0.04(syst.) 



The Belle measurement of this ratio for the cc events [ 6 1 
is consistent with the previous one: 0.55 ± 0.13(stat.) ± 
0.08(syst). 

For what concerns the D* (2458) + mass determination 
from B decays, Belle obtaines a value of 2459.2 ± 2.6 
MeV/c 2 (with combined statistical and systematic uncer- 
tainties) which is compatible with the measurements from 
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Figure 9. The Belle Collaboration I IIOI . Invariant mass distribu- 
tions for D*j candidates produced in the reaction B — » DD' s] and 
reconstructed in the three decay modes (a) D s n°, (b) D^n and (c) 
D s y. Points with errors represent the candidates associated to the 
AE signal region whereas hatched histograms those falling into 
AE sidebands; AE = (2, Ef) - Ej, eam where Ei, mm = is the 

beam energy and £, are the energies of the decay products of the 
B meson in the center-of-mass frame. 



continuum events either by CLEO and by BaBar and them- 
selves. Using B reconstruction to obtain mass values for 
both states may be the cleanest method since the feed 
across corrections are eluded. 



4 Properties of the new D sJ states and dis- 
cussion about their nature. 

There is no question about the existence of such two new 
states but their nature need still to be investigated. 

The masses of the new states appear surprisingly low for 
quark model spectroscopy; specifically D* y (2317) + is be- 
low DK threshold and D* y (2458) + falls below D*K one. 
If they are interpreted as ordinary cs states, this fact 
forces them to decay mainly via isospin-violating transi- 
tion making their widths quite narrow. Isospin is violated 
in £>; y (2317) + -> D>° and £r y (2458) + -> D* + n° since the 
and its excitations are / = 0, while the pion is / = 1. 
A possible mechanism is the decay through a virtual rj fol- 
lowed by 77 - 7T° mixing HI II . 

After the announcement of the D* y (2317) + discovery se- 
veral possible interpretations have appeared: since the 
D* y (2317) + contradicts the current models of charmed 
mesons spectroscopy either these models need modifica- 
tion or the observed state might have a non-gi/ nature. 
Among the "exotic" explanations, a DK molecule [ 1121 . 
a four-quark state [ 13 1 and a Dn atom [ 14 1 have been 
proposed. Among "ordinary" explanations there have been 
attempts to revise f 1 1 51 or modify [Q6| potential models 
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and to couple chiral perturbation theory with HQET M17I [ 
1 8 1 [ 1191 [ 1201 . Recently it has been even proposed that 
these states are a mixture of cs and four-quark states f 12 II . 

On experimental side, in order to clarify the nature of these 
states it is necessary to determine their quantum numbers 
and decay branching fractions, particularly those for radia- 
tive decays. The spin-parity of these states can be inferred 
from their decay modes. 

For a parity-conserving decay of D* y (2317) + to D^n 
only a natural spin-parity series is allowed, that is J p = 
{0 + , 1~, 2 + , ...}. J p = + is suggested by considering at the 
same time that (1) the mass is low compared to that of 
the D sl (2536) + and D* y (2573) + and (2) the efficiency cor- 
rected helicity distribution is consistent with being flat [ 
13, as expected for a spin-zero particle or for a particle 
of higher spin that is produced unpolarized. If J p = + 
the D* y (2317) + cannot decay to D + S y, whereas the decay 
to D* s + y is allowed by parity and angular momentum con- 
servations. Moreover it could not decay into three pseu- 
doscalars. The absence of the D + S y [[3) [0 [[6| and of 
the D + s ii + n~ [ 5 1 [ 6 1 decay modes makes a J p = + assign- 
ment most likely. CLEO [ 22 1 has also looked for neutrally 
charged states in D^n + and doubly charged states in 
without finding any signal. This is an argument against any 
molecular interpretation. 

On the other hand for the D* y (2458) + an unnatural spin- 
parity is suggested and thus this state is consistent with be- 
ing the missing J p — 1 + P-wave state. The observation of 
the D; y (2458) + -> D+y decay [|6| [[lO) rules out the J = 
assignment and favours a 1 + interpretation. Angular distri- 
butions from the B decays to this state indicate directly that 
spin one is strongly preferred Ml 01 . Moreover J — 1 + al- 
lows a strong, isospin-conserving but OZI-suppressed, de- 
cay to D*n + n~. Indeed this di-pion decay mode has been 
observed by Belle [ 6|. The apparent absence of the decay 
D* y (2458) + -» D* y (2317) + y [0 [|5) may indicate that the 
electromagnetic decay mechanism cannot compete with a 
likely strong but isospin-violating process resulting from 
rj - n° mixing. 

In summary it must be pointed out that most of the proper- 
ties of the new states, D* y (2317) + and ZT y (2458) + , can be 
explained if these particles are cs states. The most likely 
assignment for their spin-parity is + and 1 + . Models based 
on chiral simmetry coupled to HQET [ 1171 [ |l8 1 predict 
that these + and 1 + are the chiral partners of the 0~ and 1" 
states. This "parity doubling" prediction is experimentally 
confirmed: the chiral mass splitting between the + and 
0~ states, that is m(D* y (2317) + ) - m{D + s ), is equal within 
experimental error to that between the 1 + and 1" states, 
namely m(D* y (2458) + ) - m{D* + ). Moreover the expected 
rates for radiative decays in the same models are mostly 
consistent with experimental measurements or limits. 

Further experimental information coupled with theoretical 
ideas can definitively shed light on the nature of these new 



mesons. 
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